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ABSTRACT 

Using microsimulations of vehicular dyuamics, we 
studied the effects of several proposed iufrastructure 
developmeuts to the meau travel delay time A and mean 
speed E of vehicles passing a busy three-road fork, par¬ 
ticularly in the non-signalized roundabout junction 
of Lower Bicutan, Taguig City, Metro Manila. We 
designed and implemented multi-agent-based microsim¬ 
ulation models to mimic the autonomous driving 
behavior of heterogeneous individuals and measured the 
effect of various proposed infrastructure developments 
on A and E. Our aim is to find out the best infras¬ 
tructure development from among three choices being 
considered by the local government for the purpose of 
solving the traffic problems in the area. We created 
simulation models of the current vehicular traffic sit¬ 
uation in the area using the mean travel times r of 
statistically sampled vehicles to show that our model 
can simulate the real-world at a significance level of 
a = 0.05. Based on these models, we then simulated 
the effect of the proposed infrastructure developments 
on A and E and used these metrics as our basis of 
comparison. We found out that the proposed widening 
of one fork from two lanes to three lanes has the most 
improved metrics at the same a = 0.05 compared to 
the metrics we observed in the current situation. Under 
this infrastructure development, the A increases lin¬ 
early {B? = 0.98) at the rate of 1.03 s, while the E 
decreases linearly {E? > 0.99) at the rate of 0.14 km/h 
per percent increase in the total vehicle volume V. 


Categories and Subject Descriptors 

K.3.7.1 [Computing Methodologies]: Artificial 
Intelligence— Multi-agent systems', K.5.1.2 [Computing 
Methodologies]: Modeling and Simulation— Model 


verification and validation; L.9.5 [Applied Com¬ 
puting]: Operations Research— Transportation 

1. INTRODUCTION 

The increasing number of vehicles that use the road 
at the same time is causing congestions on road net¬ 
works especially at road sections where nominal widths 
change such as intersections, three-way forks, and 
roundabouts. One real-world example of such width 
change is the Bicutan Roundabout (BR) in Taguig 
City, Metro Manila. Commuters passing through this 
roundabout experience delay in travel everyday, espe¬ 
cially during rush hours on week days, because of 
vehicular congestion in the area. BR is located at the 
Upper Bicutan area and has an inscribed diameter of 
34 m. It is bounded to the East by the Department 
of Science and Technology (DOST) campus, to the 
North by the Philippine National Construction Corpo¬ 
ration (PNCC) campus, to the West by the Philippine 
National Railways (PNR) crossing, and to the South 
by North Daanghari. BR serves as the T-intersection 
between General Paulino Santos Avenue (GPSA) which 
lies along the ENE-WSW line and SLEX’s East Ser¬ 
vice Road (ESR) which lies along the NEN-SWS line 
of the area (Eigure[T]). The nominal center of BR is 
approximately 100 m WSW from DOST’s main gate 
along GPSA, 150 m SWS of the ESR and San Martin 
de Porres fork, and 106 m ENE of PNR crossing. 
Jeepneys, buses, trucks, taxis, AUVs, motorcycles, tri¬ 
cycles, and bicycles are the usual type of vehicles that 
pass the area. Upper Bicutan, a barangay under the 
urbanized City of Taguig, is situated on the western 
shore of Laguna de Bay and is bordered in the south 
by Muntinlupa City, in the southwest by Paranaque 
City, in the west by Cainta, Rizal, in the northeast by 
Taytay, Rizal and in the north by Makati City, Pasig 
City and the Municipality of Pateros. According to the 
2010 census, 62,570 out of the city’s total population of 
613,343 are residents of Upper Bicutan [la] . 

The thrust of the city government administration is to 
make Taguig City a “Premier City” in the Philippines. 
However, the physical and financial development of a 
city is greatly dependent, among others, on the effi¬ 
ciency of its transportation support system, particularly 
how the city manages the vehicular flow along its road 


Figure 1: The area map of Upper Bicutan, Taguig City showing the relative location of BR [^. 


networks. The efficient flow of vehicles carrying hnman 
and non-human capitals results in the efficient flow of 
resources. This translates to a productive citizenry and 
successful businesses. Because of this, the city govern¬ 
ment must make its move in pursuing efficient and prac¬ 
tical steps to at least mitigate, if not totally eradicate, 
traffic congestion which must result into improved vehic¬ 
ular flow. Vehicular flow may be quantified by the fol¬ 
lowing metrics and we say that improvement has been 
achieved if there is a statistically significant difference 
the measurement at a particular significance level a as 
compared to the current situation: 

1. Average vehicular delay time while in BR (A) - 
This metric measures the total amount of time 
that the vehicle is in the stopped position (i.e., 
when the vehicle has a speed of 0 km/h)-, and 

2. Average vehicular speed while in BR (E) - This 
metric measures the mean speed of a vehicle while 
traveling within the BR area. 

Ideally, we wanted A to approximate zero while we 
wanted E to be near but not more than the legal max¬ 
imum speed. In realistic situations, however, A > 0 and 
E is even way below the legal minimum speed. These 
two metrics complement each other because for some 


agents such as those who drive big trucks, an almost 
zero A and a very low E is much more desirable than 
a A > 0 and a relatively high E. This is because it 
requires much more fuel to move the heavy vehicles from 
a full stop than to keep them moving at low E. For 
agent drivers of light- and medium-weight vehicles, a 
A > 0 but a very high E is desirable because fuel burnt 
at low E is greater than when idle (for at least within 
some reasonable time). 

There are two ways the city can make its move: 
(1) Improve the infrastructure of the road network 
that will best optimize the A and the E, or (2) imple¬ 
ment responsive but efficient time-dependent traffic 
schemes that will significantly improve both metrics. 
Both moves are very costly if the government is to 
try all possible schemes and infrastructure just for the 
simple reason of finding the optimal one. 

In terms of finding the optimal infrastructure improve¬ 
ment for BR, one should not do it via a trial-and-error 
way because it means the following: (1) building the 
infrastructure, (2) letting it tried by the commuters for 
a reasonable amount of time just to statistically compile 
the A and the E metrics, and then (3) removing it to 
build the next one. By the time the government has 
compiled the data it needed, it has already spent sig- 



nificant amount of money, time, and other resources. 
Besides, the data collected might have already been 
biased due to what we call “data collection dependency” 
brought about by the memory of the agents in their 
experiences. This means that for a trial-and-error with 
n iterations, the data collected on the {i + l)th iterate 
might have been interactively affected by the improve¬ 
ment made during the ith iterate, unless of course, after 
every iterate we revert back to the current infrastruc¬ 
ture for the same amount of time just so both the drivers 
and the commuters will forget their ith experiences. 

In terms of implementing responsive yet efficient traffic 
schemes, most local governments have the habitual 
liking of doing it in the trial-and-error way, as expe¬ 
rienced by all of this paper’s authors. If for example 
the local Traffic Management Office (TMO) is consid¬ 
ering n traffic schemes, they will: (1) implement the 
ith scheme for a t period of time (where t is empir¬ 
ically selected to be long enough for the drivers and 
the commuters to adapt to it), (2) compile the A 
and the E metrics, and (3) repeat the process for the 
{i + l)th scheme. The local TMO will then compare 
the schemes and decide based on the metrics collected 
without considering (maybe even knowing) that there 
is an interactive effect on the {i + l)th scheme the 
experiences and adaptation of the driver and commuter 
agents during the implementation of the ith scheme. 

One efficient way that may be used to try out ideas 
and schemes in solving traffic congestions is to use 
computer-based simulation, a model created to reflect 
real-world situations. Once the model is constructed, 
solutions can be tried out to determine its impact on 
traffic congestion [^, without the need for a costly 
real-world trial-and-error implementation. In this 
study, we used a multi-agent framework to simulate the 
object-following behavior of heterogeneous agents that 
share a road network. The road network reflects the 
scaled-down dimension of a real-world network, such 
as that of the BR, servicing scaled-down dimensions 
of real-world vehicles. Our model can be character¬ 
ized as microscopic, stochastic, discrete time-step, and 
behavior-based. The model is microscopic because it 
simulates the behavior of the vehicles themselves while 
on the road in the presence of, and interacting with, 
other vehicles. The model is stochastic because it is 
based on random events that follow some distribu¬ 
tions that we observed in the real-world. The model 
is discrete time-step because at each unit time, the 
agents concurrently change their respective positions, 
directions, and speeds driven only by their respective 
simulated behaviors in response to their respective 
interactions with the current state of the road network. 
The model is behavior-based because the traffic flow 
algorithms that we used are based on a psycho-physical 
car following model for vehicles moving along the length 
of the road [i,i. We used a rule-based algorithm for 
lateral movement that simulates the lane-changing 
behaviors of the agents 0,0. 

We present in this paper the result of our experiments 
in searching for the optimal infrastructure improvement 
for BR from among three choices that are being consid¬ 


ered by the local government. We present our results in 
searching for the optimal traffic schemes for the current 
infrastructure in another paper [I^. Our study had 
three objectives: (1) To show that our model reflects 
the current real-world scenario at a statistical signifi¬ 
cance of a = 0.05 using the mean travel time r of the 
sampled vehicles as metric, (2) to use the model to sim¬ 
ulate the effect on A and on E of three infrastructure 
developments (ID) in BR currently being considered by 
the local government, and (3) simulate the effect on A 
and on E the 10%, 50%, and 100% increase in vehic¬ 
ular volume under the best ID. Based on our results, 
we found out that: (1) The average travel time of the 
simulated vehicles while at the BR is not significantly 
different at a = 0.05 from the observed average travel 
time of the sampled real-world vehicles, (2) the best 
ID is when the South bound lane of the ESR is widen 
up to three lanes continuing to the west bound lane of 
GPSA from BR to PNR, and (3) the A increases linearly 
{FI? = 0.98) at the rate of 1.03 s per percent increase 
in the total vehicle volume V, while the E decreases lin¬ 
early {I? > 0.99) at the rate of 0.14 km/h per percent 
increase in V. 

2. MODEL 

A computer model that controls an agent’s behavior 
(i.e., the driver) in terms of its reaction to the vehicle in 
front of it while at the same lane is called a car-following 
model 0- An agent-driven vehicle is said to be fol¬ 
lowing when it is constrained by the speed of the moving 
vehicle in front of it, and that driving at the agent’s 
desired speed will lead to a collision. When a driver 
agent is not constrained by another vehicle, it is said 
to be free and travels, in general, at its desired speed. 
The actions of the following agent is defined by the 
agent’s acceleration, although some models, for example 
the car-following model developed by Gipps 0], define 
the agent’s actions through the agent’s speed. Some 
car-following models only describe the agents’ behavior 
when they are currently following another vehicle, while 
other models determine the agents’ behavior in all situa¬ 
tions. We believe, however, that to model the real-world 
decision-making capabilities of human drivers, a car¬ 
following model should identify both of the following: 

1. The current state S the vehicle is in; and 

2. What actions A are desirable at this state. 

Most car-following models use several states to describe 
the following agents’ behavior. Most models use: 

1. S'/: A state for free driving, where the vehi¬ 
cles are unconstrained and the respective driver 
agents try to achieve their desired speeds (subject 
of course to pertinent legal speed limits of the road 
network); 

2. Sn'. A state for normal following, where the fol¬ 
lowing agents adjust their respective speeds with 
respect to the speeds of the vehicles in front of 
them; and 


3. Se'. A state for an emergency deceleration, 
where agents try to decelerate to avoid a collision 
with the vehicle in front. 


Throughout the paper, we used the notations summa¬ 
rized in Table[T]to describe the kinematic quantities and 
model outputs. 


2.1 Sf Free Driving State 

The agents at the Sf state try to accelerate or decel¬ 
erate to achieve its current desired speed. If the current 
speed Vi of the ith vehicle is higher than its desired 
speed the agents uses the normal deceleration 

rate (—to slow down to the desired speed. If 
Vi < , the agents use its maximum acceleration 

rate to reach at the shortest time possible. 

The — and are parameters of the simula¬ 

tion and they are dependent on the type of vehicle the 
agent is driving. The acceleration rate of the ith agent 
is expressed as 


ai = 


The time ti it will 
desired speed is 


ti = 


max 

Vi < 

0 

Vi = uf 

norma.1 

(Ij^ 

Vi > uf 

take for 

the ith vehicle to achieve its 

Vi 

^ max 

Vi < 


Vi > wf 

Q^normal 


3.1 Interview with Stakeholders 

We conducted interviews in order for us to understand 
the current plans of the two local governments that have 
jurisdictions over the BR: the Taguig City Traffic Man¬ 
agement Office (TMO) and the Metro Manila Develop¬ 
ment Authority (MMDA). The interview served as some 
sort of leveling off with the intended stakeholders. We 
also wanted the result of this research endeavour to be 
put to use, particularly because two of the three authors 
of this paper pass this area on a daily basis. Thus, 
if either or both local governments will implement the 
result of our research, we will greatly benefit from such. 
We also asked permission from the two governments to 
allow us to conduct scientific observations so that we 
can quantify the current vehicular flow in BR. 


3.2 Observation of Current Situation 

We employed the help of twelve enumerators whom we 
assigned to different identified points within the BR 
junction. Each point of entrance to and exit from BR 
has designated persons to record the time of entrance 
and exit of the vehicle, and the type of vehicle that 
passed through. Each vehicle is identified by their tag 
or plate number. We recorded a series of actual obser¬ 
vations during peak and non-peak periods for us to be 
able to get a statistically accurate data and determine 
the highest and lowest traffic volume. From these obser¬ 
vations, we were able to obtain the mean travel time r of 
sampled vehicles, as well as the respective distributions 
of each vehicle type. To identify the specific points in 
the BR, we divided the area into six routes as follows: 


2.2 Sn Normal Following State 

While at the normal following state 5'„, the accelera¬ 
tion rate of the ith vehicle is given by an a^mmetrical 
Gazis-Herman-Rothery (GHR) model [3, 0, [3- The 
acceleration is computed as 


where r^, s^, and are model parameters. The 
parameters r^, s"*", and t"*" are used if Vi < Vi-i, while 
the parameters r~, s~ , and t~ are used if Vi > Vi-i. 
Notice here that the final acceleration of the ith vehicle 
is given as max{ai,ai_i}. 

2.3 Emergency State 

Under this state, the agents use a deceleration rate that 
prevents collision and extends Sx. This deceleration rate 
is given by 


1. Route 1 is the route from PNR to DOST Campus 
Eastbound along the GPSA passing through the 
BR; 

2. Route 2 is the route from PNR to PNCC 
Campus, Eastbound along GPSA from the PNR 
to BR, and then turning left through the BR, and 
up to the Northbound lane of ESR going towards 
the PNCC Campus; 

3. Route 3 is the route from DOST campus to PNR 
Westbound along the GPSA passing through the 
BR; 

4. Route 4 is the route from the DOST Campus to 
the PNR Crossing, Westbound along the GPSA 
up to the BR, and then turning right through the 
BR, and up to the Northbound lane of ESR going 
towards the PNCC Campus; 


• r normal — 4 

minjai ,ni-i - Vi > Vi-i 

min{ar™',ai-i + 0.250?°™"' * Vi < Vi-i 

3. METHODOLOGY 

We present in this section the activities that we per¬ 
formed in this study. We also present here the metrics 
that we measured, as well as the statistical analysis that 
we employed to analyze the results of the study. 


5. Route 5 is the route from PNCC Campus to 
PNR, Southbound along the ESR, turning right 
through BR, and then Westbound along the 
GPSA to the PNR; and 

6. Route 6 is the route from PNCC Campus to 
the DOST Campus, Southbound along the ESR, 
turning left through BR, and then Eastbound 
along the GPSA to the DOST Campus. 





Table 1: Alphabetical listing of notations, mathematical variables, and abbreviations used in this 
paper including their respective descriptions. 

Notation Description 

5x Space headway between x„-i and Xn. in m 
5v difference in speed between Xn-i and Xn in m/s 
Un Acceleration of the nth vehicle in m/s^ 

Ln length of the nth vehicle in m 

s„ Effective length of the nth vehicle in m 

T Reaction time in s 

Vn Speed of the nth vehicle in m/s 

^desired Deshed Speed of the nth vehicle in m/s 
Xn Longitndinal position of the nth vehicle in m 
a Statistical level of significance 

A Mean delay time per vehicle in s 

E Mean vehicular speed in m/s 

T Mean vehicular travel time s 

To T of the sampled vehicles from observation 
Ts Mean r of the simulated vehicles 

V Total volume of vehicles 

ANOVA Analysis of Variance using the F-Statistics 
BR Bicutan Roundabout 

DOST Department of Science and Technology 

ESR East Service Road 

GPSA General Paulino Santos Avenue 

ID Infrastructure Development 

MMDA Metro Manila Development Authority 

PNGC Philippine National Construction Corporation 

PNR Philippine National Railways 

PUP-T Polytechnic University of the Philippines-Taguig 
SLEX Southern Luzon Expressway 

TUP-T Technological University of the Philippines-Taguig 
TMO Traffic Management Office 


3.3 Simulation of Current Situation 

We conducted a replicated microsimulation study of the 
vehicular flow under the current BR using the data on 
respective distribution of vehicles by type. We repli¬ 
cated the study to n = 10 and computed the mean r. 
We wanted to know if the model can statistically repro¬ 
duce the observed vehicular traffic data. Statistically, 
the respective differences of the r between the observed 
ijo) and the mean simulated (ts) runs must not be dif¬ 
ferent from zero at a = 0.05. We used analysis of vari¬ 
ance statistics (ANOVA) to evaluate two hypotheses, 
the null hypothesis Hq and the alternative hypothesis 
Ha as follows: 

Hq-. There is no significant difference between To and Ts 
at a = 0.05. 

Ha- There is significant difference between To and Ts at 
the same a. 

3.4 Microsimulations of Proposed Devel¬ 
opments 

We created the respective replicated n = 10 microsim¬ 
ulation studies of the vehicular flow using the observed 
current vehicular distribution when under the different 
planned infrastructure developments (ID) of the BR. 
These IDs are based on the G6 Project blueprint, the 
medium-term development plan by the Gity Mayor of 


Taguig Lani Gayetano, and DOST’s Monorail Project 
which is currently underway. The proposed develop¬ 
ments are: 

1. Infrastructure Development 1 (IDl): Widening of 
the GPSA into three lanes but only those lanes 
towards the BR. That is the east bound lane from 
PNR to BR and the west bound lane from DOST 
to BR. 

2. Infrastructure Development 2 (ID2): Widening 
of the GPSA into three lanes but only the east- 
bound lane from PNR to BR continuing from BR 
to DOST gate. 

3. Infrastructure Development 3 (ID3): Widening of 
the south bound lane of the ESR into three lanes 
up to the BR continuing to the west bound lane 
of GPSA from BR to PNR. 

We computed the respective means of A and S for all 
infrastructure improvements, including that of the cur¬ 
rent infrastructure we termed here as IDO (i.e., from 
IDO through ID3) and conducted an analysis of vari¬ 
ance (ANOVA) statistics to find out whether the respec¬ 
tive means are significantly different from each other at 
a = 0.05. That is, we want to test two hypotheses, the 
null hypothesis Hq and the alternative hypothesis Ha as 
follows: 



Hq-. The absolute pairwise differences |Aido — Aidi| = 
• • • = |Aid1 — A1D2I = |Aid1 — AidsI = |Aid2 — AidsI = 
0, where the = sign in this sentence wonld mean “not 
signihcantly different from” or “statistically equal to.” 

Ha’- Any of the following respective pairwise absolnte 
differences is true: |Aido — Aidi| > 0, ... , |Aidi — 
A 1 D 2 I > 0, |Aidi — AidsI > 0, or |Aid 2 — Aids! > 0- 

We have the same set of hypotheses for E. 

3.5 Increased Traffic Volume at Best ID 

We rerun the respective replicated microsimulation 
studies under fDO, IDl, ID2, and IDS but with respec¬ 
tive increases of 10%, 50%, and 100% in V, making sure 
that the distribution by vehicular type is preserved. 
Here, we wanted to know if the benefits of the proposed 
improvements will be carried with the increase in V and 
if so, up to how much increase. The assumed increase 
in V is just a natural reaction of the driving agents 
when there is a perceive improvement in the current 
situation. The improvement of the vehicular flow due to 
infrastructure improvements will attract more agents, 
increasing V in the area, and thereby might degrade 
the expected designed benefits of the improvement in 
the long run. 

4. RESULTS AND DISCUSSION 

4.1 Actual Current Vehicular Flow 

Tables [2] and [3] summarize the distribution of vehicle 
per route and per vehicle type, respectively, that we 
compiled during the actual observation. Route 4 has 
the most number of total vehicles that passed during 
the observation period with 23.79%. This is followed 
by Routes 5, 1, and 2, in non-increasing order. Route 
4, undeniably has the most number of vehicles since 
Route 4 is the only route that may be used by residents 
from the Lower Bicutan area who are trying to go to 
the Metropolitan Manila Area via the C-5 Diversion. 
Aside from the Lower Bicutan residents, the route is 
also extensively used by students and employees of sev¬ 
eral Higher Education Institutions (HEIs) and govern¬ 
ment agencies such as the Technological University of 
the Philippines-Taguig (TUP-T), the Polytechnic Uni¬ 
versity of the Philippines-Taguig (PUP-T), the LGU- 
funded Taguig City University, and DOST’s national 
offices and laboratories. Route 5 is extensively used 
by trucking companies with medium-length trucks car¬ 
rying heavy cargoes. These vehicles are coming off from 
the C-5 Diversion going to the Southern Luzon Area 
via the SLEX. Route 1 is obviously used by commuters 
coming from the SLEX area going to the residential 
Lower Bicutan or to any of the HEIs and government 
agencies. Route 2 is the reverse of Route 5, where 
mostly medium-length trucks pass through carrying car¬ 
goes from the Southern Luzon area via SLEX going to 
Metropolitan Manila Area through the C-5 Diversion. 
From Table [31 motorcycles and cars combined account 
for about 70% of the vehicles that use BR. This is fol¬ 
lowed by bicycles, jeepneys and vans. The heavier vehi¬ 
cles, such as the many-wheelers and buses only account 
for about 4% combined. This distribution by type is 
actually expected because the area is used extensively 


by residents and commuters rather than by non-human 
cargoes. Table U shows the mean vehicular dimensions 
by vehicle type. We used these dimensions to scale the 
respective agents in our simulations. 

Table 2: Distribution of the total number of 
vehicles that passed through the BR during the 
observation period. 

Route Number of Vehicles 

Percentage (%) 

Route 1 

778 

20 

Route 2 

719 

18 

Route 3 

524 

13 

Route 4 

934 

24 

Route 5 

815 

21 

Route 6 

156 

4 

Total 

3,926 


Table 3: Distribution of the vehicles per type. 

Type 

Percentage (%) 

Motorcycle 


38.30 

4x8 wheeler 


0.64 

4x6 wheeler 


2.75 

Van 


5.88 

Jeepney 


9.60 

Car 


31.73 

Bus 


0.38 

Bicycle 


10.72 

Total 


100.00 

Table 4: The average 
vehicles per type. 

dimension of observed 

Vehicle 

Length 

Width 

Type 

(m) 

(m) 

Motorcycle 

2.00 

1.5 

4x8 wheeler 

6.59 

1.5 

4x6 wheeler 

5.41 

1.5 

Van 

5.50 

1.5 

Jeepney 

4.00 

1.5 

Car 

4.50 

1.5 

Bus 

11.54 

2.5 

Bicycle 

1.45 

0.5 


4.2 Evaluation of the Model against the 
Observed 

Table \B summarizes the statistics we computed after 
conducting ANOVA on Tq and Ts- Since we see that the 
F statistics ap = 0.6873 > a = 0.05, then we accept 
the null hypothesis Ho and we say that the difference 
between the mean sampled observed data To and the 
mean simulated data Ts is not significantly different from 
zero at a « 0.05. Thus, our microsimulation model was 
able to mimic the driver’s behavior in the real world 
with a guarantee of being correct 1 — a = 0.95 of the 
time. 

4.3 Comparison of the IDs 

Tables [5] and [7| show the respective ANOVA tables of 
the mean A and mean E for IDO, IDl, ID2, and ID3 
where we find that Op < a and Op < a. Based on the 



















Table 5: The ANOVA table comparing To and Ts using the F statistics. SOV means Source of Variation 
and DF means Degree of Freedom 


SOV 

DF 

Sum of Squares 

Mean Square 

F 

Qip 

Replication 

9 

1,633.05 

181.45 

39.21 

< 0.0001 

To VS. Ts 

1 

0.80 

0.80 

0.17 

0.6873 

Error 

9 

41.65 

4.63 



Total 

19 

1,675.50 





computed Op and Op, we accept the alternate hypoth¬ 
esis Ha saying with confidence 1 — ap = 0.9999 that 
at least one of the absolute pairwise difference between 
AiDa: and Aid;/, and between SiDa: and Sidh, Vx 7 ^ y, is 
significantly greater than zero. Figures [5] and [3] show the 
respective mean A and mean E of IDs;, Va; = 0,..., 3. 
Based on the DMRT, we see that IDl will not improve A 
nor E. In fact, both A and E under IDl is no different, 
statistically, than the A and E of the current infrastruc¬ 
ture, and thus, implementing IDl will just be a waste of 
taxpayer’s money. We also see that ID2 has significantly 
improved A and E compared to the current infrastruc¬ 
ture. However, IDS has the most significant improve¬ 
ment in terms of both A and E compared to the current 
infrastructure. From Figure [3] we see that the respec¬ 
tive E of ID2 and IDS are not statistically different from 
each other, so either ID can provide the best E. Based 
on the two metrics, however, we say that IDS can pro¬ 
vide the improvement that we desire for the vehicular 
flow. 


IDO IDl ID2 ID3 

Infrastructure Development 

Figure 2: The respective mean A of the var¬ 
ious IDs with DMRT grouping. Means with the 
same letter are not statistically different from 
each other by DMRT. 

4.4 Effect of Increased Traffic Volume 

Figures|4]and[5]show the respective effect of increasing V 
from 10% to 100% on A and E while under IDS. As V 
is increased one percent at a time, A increases linearly 
with = 0.98 at a rate of 1.03 s while E decreases 
linearly with H? > 0.99 at a rate of 0.15 km/h. At this 
linear rate of decrease, all vehicles will come to a halt 
(i.e., E = 0) at V « 180%. 

5. CONCLUSION 

In this paper, we presented a multi-agent-based 
microsimulation of a real-world vehicular traffic based 
on the existing psycho-physical car following model 



Infrastructure Development 

Figure 3: The respective mean E of the var¬ 
ious IDs with DMRT grouping. Means with the 
same letter are not statistically different from 
each other by DMRT. 

for simulating the movement of vehicles along the 
lane [ 3 , @| and a rule-based algorithm for simulating 
the lane-changing behaviors of the vehicles 0 , 0 . 
We simulated a real-world infrastructure involving six 
routes merging into a three-road fork via an unsignal¬ 
ized roundabout. This roundabout is the BR located in 
Taguig City, Metro Manila. From the observed distri¬ 
bution of vehicles per route and per vehicular type, we 
simulated the vehicular traffic of the current situation 
and collected the r metric. Based on the ANOVA, we 
found out that To and Ts are not statistically different 
from each other based on the F statistics. This means 
that our simulation model can represent the real-world. 
We then simulated three proposed infrastructure devel¬ 
opments (IDl, ID2, and IDS) and found out that IDS 
offers the least A and at the same time the fastest E. 
We increased the V from 10% to 100% under IDS and 
found out that A increases linearly with V, while E 
decreases linearly. Based on these results, we provide 
the following conclusions: 

1. Our microsimulation model can stastically repre¬ 
sent the real-world; 

2. The best ID for BR is IDS; and 

3. Under IDS, V has a positive linear effect on A and 
a negative linear effect on E. 

6. RECOMMENDATIONS 

We recommend to the stakeholders within BR, such 
as the City Government of Taguig through its TMO, 
as well as the Engineering Section of the MMDA, the 
residents of Lower Bicutan, the public and private 
















Table 6 : The ANOVA table comparing the different IDs in terms of A using the F statistics. SOV 
means Source of Variation and DF means Degree of Freedom 

SOV DF Sum of Squares Mean Square F Op 

Replication 9 2,328.63 258.74 5l2 0.0004 

ID 3 3,556.14 1,185.38 23.47 < 0.0001 

Error 27 1,363.39 50.50 

Total 39 7,248.15 


Table 7: The ANOVA table comparing the different IDs in terms of E using the F statistics. SOV 
means Source of Variation and DF means Degree of Freedom 


SOV DF Sum of Squares Mean Square F ^ 


Replication 

9 

14.73 

1.64 

5.63 

0.0002 

ID 

3 

11.83 

3.94 

13.56 

< 0.0001 

Error 

27 

7.85 

0.29 



Total 

39 

34.41 






Figure 4: The respective mean A of IDS when V 
was increased from 10% to 100%. The line rep¬ 
resents the regression with A = 1.03 x V + 72.64 
= 0.98). 

employees within BR such as the DOST, the TUP-T, 
and the PUP-T, and the private sector, particularly 
SM Corporation, to help call for the implementation 
of infrastructure developments described in ID3. We 
caution, however, that we may be able to find a better 
vehicular flow if hybrid improvements of IDl through 
ID3 are considered. We further caution that it may 
be cheaper to implement signalized traffic schemes 
together with the ID. As we have mentioned before, 
we have already reported our findings with various 
traffic schemes elsewhere [l^, while we are currently 
working @] on the interactive effects of ID and traffic 
schemes to both A and E. 

7. FUTURE WORK 

Several scientific research efforts have already forked out 
from this initial investigation. These are: 

1. Studying the effects of installing railings in busy 
sidewalks; 

2. Pedestrian dynamics for a very, very large aca¬ 
demic building beyond its designed full capacity; 



Figure 5: The respective mean E of IDS when V 
was increased from 10% to 100%. The line rep¬ 
resents the regression with E = —0.15 x V -I- 26.93 
{R^ > 0.99). 


3. Evaluation of evacuation plans for Tsunami- 
vulnerable shoreline communities; 

4. Designing re-routable and adaptive traffic plans 
for time-responsive traffic schemes; and 

5. Incorporating the “greedy” behavior of Jeepney 
and tricycle drivers in the car-following model. 
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